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Abstract 

Using the spectral distribution associated with the adjacency matrix of graphs, we 
introduce a new method of calculation of amplitudes of continuous-time quantum walk 
on some rather important graphs, such as line, cycle graph C^, complete graph Kn, graph 
Gn, finite path and some other finite and infinite graphs, where all are connected with 
orthogonal polynomials such as Hermite, Laguerre, Tchebichef and some other orthog- 
onal polynomials. It is shown that using the spectral distribution, one can obtain the 
infinite time asymptotic behavior of amplitudes simply by using the method of stationary 
phase approximation (WKB approximation), where as an example, the method is applied 
to star, two-dimensional comb lattices, infinite Hermite and Laguerre graphs. Also by 
using the Gauss quadrature formula one can approximate infinite graphs with finite ones 
and vice versa, in order to derive large time asymptotic behavior by WKB method. 
Likewise, using this method, some new graphs are introduced, where their amplitude 
are proportional to product of amplitudes of some elementary graphs, even though the 
graphs themselves are not the same as Cartesian product of their elementary graphs. 
Finally, via calculating mean end to end distance of some infinite graphs at large enough 
times, it is shown that continuous time quantum walk at different infinite graphs belong 
to different universality classes which are also different than those of the corresponding 
classical ones. 

Keywords: Quantum walk, Continuous-time quantum walk, Spectral dis- 
tribution. Graph. 
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The theory of Markov chains and random walks on graphs is fundamental to mathematics, 
physics, and computer science [1, 2, 3] as it provides a beautiful mathematical framework to 
study the stochastic process and its applications. Among the known examples of the applica- 
tions we can mention Monte Carlo methods in statistics, the theory of diffusion in statistical 
physics, and algorithmic techniques for sampling and random generation of combinatorial 
structures in computer science (based on rapid mixing of certain Markov chains). Two perva- 
sive algorithmic ideas in quantum computation are Quantum Fourier Transform (QFT) and 
amplitude amplification (see [4]). Most subsequent progress in quantum computing owed much 
to these two beautiful ideas. But there are many problems whose characteristics match neither 
the QFT nor the amplitude amplification mold (e.g., the graph isomorphism problem), where 
this begs for new additional tools to be discovered. 

A natural way to discover new quantum algorithmic ideas is to adapt a classical one to the 
quantum model. An appealing well-studied classical idea in statistics and computer science is 
the method of random walks [5]. Recently, the quantum analogue of classical random walks 
has been studied in a flurry of works [6, 7, 8, 9, 10, 11]. The works of Moore and Russell [10] 
and Kempe [11] showed faster bounds on instantaneous mixing and hitting times for discrete 
and continuous quantum walks on a hypercube (compared to the classical walk). 

A study of quantum walks on simple lattice is well known in physics(see [12]). Recent 
studies of quantum walks on more general graphs were described in [6, 7, 9, 13, 14]. Some of 
these works studies the problem in the important context of algorithmic problems on graphs 
and suggests that quantum walks is a promising algorithmic technique for designing future 
quantum algorithms. 

Several important classes of graphs studied in classical random walks include the binary 
n-cube, the circulant graphs, and the group-theoretic Cayley graphs. The binary n-cube and 
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circulant graphs are important in the study of interconnection networks and complexity of 
Boolean function, and Cayley graphs capture strong gorup-theortic ingredients of important 
problems, such as graph isomorphism. Since most of these graphs are regular, the classical 
random walks on them are known to converge or to mix towards the uniform stationary 
distribution. The mixing properties of continuous-time quantum walks on the same graphs 
were found to exhibit non-classical behavior [10, 15, 16, 17]. 

Here in this work we have used spectral distribution associated with the adjacency matrix, 
of some particular graphs which possess quantum distribution ( for more details see Ref.[19]) 
in studying continuous-time quantum walk. Spectral distribution helps us to give a general 
expression for the amplitude of observing the continuous-time quantum walker on a given 
site in terms of the integral over product of the spectral distribution and some polynomials 
(in most cases over the well-known orthogonal polynomial; for the details see the appendix 
II). We use isometry between the orthogonal polynomials (Hermite, Laguerre and Tchebichef) 
and interacting Fock space associated with the infinite graphs, representing raising (lowering) 
part adjacency matrix, with first order differential operator, and obtain the amplitude prob- 
ability for infinite graphs associated to orthogonal polynomials. In obtaining the amplitude 
probability for many of the finite (complete Kn, complete cycle, G„, finite path)or infinite 
graphs (line, associated graphs of orthogonal polynomials Hermite and Laguerre) we also study 
the behavior of graph for large vertices, and the asymptotic behavior of the walks. By using 
direct product graphs, we obtain the new graphs and define an approach for calculating the 
amplitude probability from the subgraphs. Also by using the Gauss quadrature formula one 
can approximate infinite graphs with finite ones and vice versa, which leads to the derivation 
of large time asymptotic hehayior continuous-time quantum walk simply by using the method 
of stationary phase. Also, we introduce graphs which have different structure, but they have 
uniform amplitude for observing particle at every level. At the end , via calculating average 
visiting strata ( sites for infinite line) which is equivalent to mean end to end distance of some 
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infinite graphs at large enough times, it is shown that continuous time quantum walk on infi- 
nite Hermite and line graphs belong to the same universality class which is different from that 
of infinite Laguerre one. 

The organization of this paper is as follows. In section 2, we give a brief outline of graphs 
and their adjacency matrices. In Section 3, we review the quantum decomposition for adjacency 
matrix of graphs, and the method for obtaining vacuum spectral distribution fi, and give some 
examples of isometry from orthogonal polynomials into interacting Fock space. Section 4 is 
devoted to the method of computing amplitude probability for continuous-time quantum walk, 
through spectral distribution of the adjacency matrix A. Section 5 is concerned with direct 
product of quantum decompositions for finite graphs. In section 6, we calculate the amplitude 
probability for continuous-time quantum walk on bunches of finite and infinite graphs. In 
section 7, in first subsection we evaluate the average visiting strata of some infinite graphs for 
large enough time, and in remaining subsections, large time asymptotic behavior of quantum 
walks on different finite and infinite graphs of section 6 are obtained by using the method of 
stationary phase(WKB) and Gauss quadrature rule. Paper is ended with a brief conclusion 
and two appendices / and //, where the first appendix consists of the proof of lemma regarding 
the equality of the amplitudes associated with vertices belonging to the same stratum and the 
second appendix contains the list of some of the orthogonal polynomials connected with some 
particular infinite and finite graphs, respectively. 



2 Graphs and its adjacency matrix 

In this section we give a brief outline of some of the main features of graphs and walk on them, 
such as adjacency matrix, stratification and orthonormal basis of strata. 
A graph is a pair G — {V,E), where F is a non-empty set and is a subset of j £ 

7^ j}- Elements of V and of E are called vertices and edges, respectively. Two vertices 
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i;j e V are called adjacent if £ E, and in that case we write i j. For a graph 

G = (y, £■) we define the adjacency matrix A — {Aij)ij^v by 

1 if i ~ J 
otherwise. 

I. 

Obviously, (i) A is symmetric; (ii) an element of A takes a value in {0, 1}; (iii) a diagonal 
element of A vanishes. Conversely, for a non-empty set V ^ a graph structure is uniquely 
determined by such a matrix indexed by V . The degree or valency of a vertex i e V is defined 

by 

k{i) = \{j e V;i^j}\, 

where | . | denotes the cardinality. A finite sequence io; ii; in £ ^ is called a walk of length n 
(or of n steps) if ik-i ~ ik for all k — 1,2, n. In a walk some vertices may occur repeatedly. 
Unless otherwise stated, we always assume that a graph under discussion satisfies: 

(a) (connectedness) any pair of distinct vertices are connected by a walk; 

(b) (local boundedness) < oo for all i E V; In fact, the examples in this paper satisfy the 
following condition which is stronger than (b): 

{b') (uniform boundedness) sitpjgy«;(i) < oo. 

Let 1^{V) denote the Hilbert space of C- valued square-summable functions on V, and G 
V} becomes a complete orthonormal basis of l'^{V). The adjacency matrix is considered as an 
operator acting in P{V) in such a way that 

A\i) = J2\j), i^V. 

Then, A becomes a self-adjoint operator equipped with a natural domain. As is easily checked, 
{h') is a necessary and sufficient condition for ^4 to be a bounded operator on ^^{V). 
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2.1 Stratification 

For i ^ j let d{i,j) be the length of the shortest walk connecting i and j. By definition 
d{i,j) — for all i & V. The graph becomes a metric space with the distance function d. Note 
that d{i,j) = 1 if and only if i ~ j. We fix a point o e F as an origin of the graph. Then, the 
graph is stratified into a disjoint union of strata: 

oo 

V^[jVk, Vk = {ieV;d{o,i) = k}. (2-1) 

k=G 

With each stratum we associate a unit vector in t^iV) defined by 

\4>k) = ^Y.\k.^)^ (2-2) 



where, i) denotes the eigenket of ith vertex at the stratum k. The closed subspace of P{V) 
spanned by {|0fc)} is denoted by r(G'). Since {|0fe)} becomes a complete orthonormal basis of 
r(G), we often write 

r(G') = E®C|0ik). (2-3) 

k 

3 Quantum decomposition 

In this section, first we review the quantum decomposition for adjacency matrix of some partic- 
ular graphs called QD graphs, and the method for obtaining the vacuum spectral distribution 
II (for more detail see Ref. [18]), then we give some examples of isometry from orthogonal 
polynomials into interacting Fock space. 

Let A be the adjacency matrix of a graph G = {V,E). According to the stratification (2-1), 
we define three matrices A'^ , A~ and A^ as follows: for i e we set 



(A+) 



Aji if j e Vk+1 
otherwise. 
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Aji if j e Vk-i 

otherwise. 

Aji if j e Vk 

otherwise. 



Or equivalently, for \k,i), 



A+\k,i)= A-\k,i)^ Yl A'\k,i)^Y\k,j), (3-4) 

for j ~ i. Since i & Vk and i ~ j then j e Vfe_i U Vit U Vit+i, where we tacitly understand that 
y_i = 0, Now with the help of [19] we define 



A^A+ + A- + A^. 



(3-5) 



This is called quantum decomposition of A associated with the stratification (2-1). Note also 
that 

(A+y^A-, {A^y^A'', (3-6) 

which can be verified easily. The vector state corresponding to |o) = \4>o), with a & V as 
the fixed origin, is analogous to the vacuum state in Fock space. According to Ref.[19], the 
< A"^ > coincides with the number of m-step walks starting and terminating at o, also, by 
lemma 2.2, [19] if T{G) is invariant under the quantum components A^, e e — , 0}, then 
there exist two Szego- Jacobi sequences {uJk}kLi and {ak}'k=i derived from A, such that 



A^\(t)k) ^ ^JuJ^\(t)k+l) , k>0 

A-|0o) = O, A-\(f)k) ^V^\(f)k-i), k>l 
= Q;fe+i|0fe), A;>0, 

where ^a}^ = 'y^+iff /^-O'), = \{i e Vfe;i j}\ for j e Vk+i and ak+i 
i^ou) = |{i e T4; i ~ j}\ for j e Vk. 



(3-7) 

(3-8) 
(3-9) 

Ko(^j), such that 
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Obviously from the relation between ujkiOtki H-{y) and | | it follows that | | — 1 > ctfe+i 
and \^''\'^^+'^ must be integer. 

Prom now on we call the graphs with the above property as kids of graphs which possess 
quantum distribution or QD type. 



4 Spectral distribution /i of the adjacency matrix A 

The spectral properties of the adjacency matrix of a graph play an important role in many 
branches of mathematics and physics . The spectral distribution can be generalized in various 
ways. In this work, following Ref.[19], we consider the spectral distribution ii of the adjacency 
matrix A: 

<A'^>= I x'^iiidx), m = 0,1,2,... (4-10) 

JR 

where < . > is the mean value with respect to a state |0o)(the ket corresponds to the ground 
stratum). By the local boundedness condition (b) of section 3 the moment sequence {< A^ > 
}m=o is well-defined[19]. Then the existence of a spectral distribution satisfying (4-10) is a 
consequence of Hamburgers theorem, see e.g., Shohat and Tamarkin [[21], Theorem 1.2]. 
We may apply the canonical isomorphism from the interacting Fock space onto the closed linear 
span of the orthogonal polynomials determined by the Szego-Jacobi sequences {{cuk} , {o.k}) ■ 
More precisely, the spectral distribution /i under question is characterized by the property of 
orthogonalizing the polynomials {Qn} defined recurrently by 

Qo{x)^l, Qi{x) ^ X - ai, 

xQn{x) = Qn+\{x) + an+lQn{x) +UJnQn-\{x), (4-11) 

for n > 1. If such a spectral distribution is unique (e.g., if the uniform boundedness condition 
ih') is fulfilled), the spectral distribution fi is determined by the identity: 

G{x)=l ^ ^gj^^^_4^ 

" jRX-y X-O^l- a- -"2 Qn{x) {r{ X - Xi' 
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where Gn{x) is called the Stieltjes transform and Ai is the coefficient in the Gauss quadrature 
formula corresponding to the roots xi of polynomial Qn{x) and where polynomials {Q!}^^} are 
defined recurrently as 

Q^i\x) ^ X - a2, 

for n > 1. 

Now if Gn{x) is known, then the spectral distribution // can be recovered from by 
means of the Stieltjes inversion formula: 

1 [V 

jji{y) — jji{x) = lim / Im{G Ju -\- iv)}du. (4-13) 

TT V — >o+ Jx 

Substituting the right hand side of (4-12) in (4-13), the spectral distribution can be determined 
in terms of x;, Z = 1, 2, the roots of the polynomial Qnix), and Guass quadrature constant 
Ai, 1 — 1,2, ... as 

I 

( for more details see Ref. [18, 20, 21, 22].) 

In the following, we show the above mentioned isometry U from the orthogonal polynomial 
into interacting Fock space in infinite graphs associated with Hermite, Laguerre and Tchebichef 
polynomials, simply by replacing the raising (lowering) part of adjacency matrix, i.e., A"*" {A~) 
by the corresponding ffist order differential operators S_(n)(>l_(n)) introduced by one of the 
authors in Ref. [23] with the following recursion relations 

B^{n)Qn-i{x) = Qn{x); A_{n)Qn{x) = E{n)Qn-i{x), (4-15) 

where E{n) — n, for Hermite, E{n) = n{n — 7), 7 > —1, for Laguerre and E{n) = v? for 
Tchebichef. Now by assuming the canonical isometry UQn{x) = -\/^|0n) between the corre- 
sponding Hilbert spaces, ^4+ = UB^{n)U* , A" — UA^{n)U* and the corresponding operators 
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for Laguerre polynomial Qn{x) = Ln{x) 

A+ ^Ui-x^+x-nW"; A' ^U(x^ -n)U\ 
ax ax 

for Hermite polynomial Qn{x) = Hn{x) 

A+ = U(-4- + x)U*; A- = U{^)U\ 
ax ax 

for Tchebichef polynomial Qn{x) = T„(a;) 

A+ = U(-(l -x^)-^ + nx)W- A- = Uiil - x')-^ - nx)U\ 
ax ax 

Finally, using the quantum decomposition relations ( 3-7,8,9) and the recursion relation 
(4-11) of polynomial (5„(x), the other matrix elements {(j)k\A^ \ 4>q) can be written as 



\A'' 



bo) = , ^ [ x"'Qk{x)ii{dx), m = 0, 1, 2, .... (4-16) 

Our main goal in this paper is the evaluation of amplitude probability for continuous-time 
quantum walk by using Eq.(4-16) such that we have 

qk{t) = (0fe|e-^^* I 0o) = , ^ / e-'^'Qk{xMdx), (4-17) 

where qkit) is the amplitude of observing the particle at level k at time t. The conservation of 
probability I]fe=o I Q'a:(^) 1^ = 1 follows immediately from Eq.(4-17) by using the completeness 
relation of orthogonal polynomials Qn{x). Also by using Eq.(4-17) and (3-7,3-8,3-9), we obtain 
the following recurrence relation for the amplitude Pfe(t): 

= V^ik+l?ik+l(^) + Q;it+i?it(^) + y/u}^qk-i{t). (4-18) 

Obviously evaluation of qk{t) leads to the determination of the amplitudes at sites belonging 
to the stratum T4, as it is proved in the appendix I, the walker has the same amplitude at 
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sites belonging to the same stratum, i.e., we have qik{t) = = 0, 1, | Vk |, where qik{t) 

denotes the amphtude walker at ith site of kth stratum. 
Obviously for finite graphs, the formula (4-17) yields 

qk{t) = , ^ E Aie-'^''Qk{xi), (4-19) 

where by straightforward calculation one can evaluate the average probability for the finite 
graphs as 

P{k) = hm 1 r I qk{t) r dt = J:AiQU^i)- (4-20) 

T-^oo I Jo UJ1UJ2 ■ ■ - UJk I 

Finally using Gauss quadrature formula [20] 

/n 
f{xMdx) = Y.Akf{xi), (4-21) 
1=1 

where the constants Ai and xi are the same as those appearing in formula(4-12), we can 
approximate an infinite QD graph corresponding to the sequence of orthogonal polynomials 
{Qo,Qi,Q2, ■■■}, with the finite n-stratum QD graph of spectral distribution obtained from 
Stieltjes transform (4-12) with constant Ak and Xk, and vice versa. Therfore, using Gauss 
quadrature formula, we can interchange formulas (4-17) and (4-19) of finite and infinite graphs. 
This approximation can be very useful in obtaining the large time asymptotic behavior of 
quantum walk on graphs with finitely many vertices. 

Since, studying the large time behavior of quantum walks on finite or infinite graphs by 
this method, naturally leads us to consider the behavior of integrals of the form 

I{t)^ f\-'*'^^''^g{x)dx a<x<b (4-22) 

Ja 

as t tends to infinity. There is a well-developed theory of the asymptotic expansion of integrals 
known as method of stationary phase approximation or WKB [24, 25], which allows us to 
determine, very precisely, the leading terms in the expansion of the integral in terms of simple 
functions of t (such as inverse powers of t) . In section 7 we will study asymptotic behavior of 
some finite and infinite graphs at distant time by using Gauss quadrature formula and WKB 
approximation. 
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5 The direct product of QD finite graphs 

In this section, we study continuous-time quantum walk on direct product of QD finite graphs. 
Let Gi,i — l,2,..,n be graphs of finite vertices with the corresponding adjacency matrices 
Ai,i — 1, 2, n. Then their direct product 

Gi ® ■ ■ ■ ® G„, (5-23) 

is a graph with the following adjacency matrix A: 

n 

A = ^/ ®- -- ®Aj®---®/ (5-24) 

where the j — th term in the sum has Aj appearing in the j — th place in the tensor product. 
In this case, for the amplitude of the ground state |0o) = \4>o^)\4''o ^) " ' \ 4>o^^)-> similar to 
Eq.(4-10), we have 

(^ol^"" \ (l)o) ^ J ■■■ J(xi+X2 + --- + XnTlJi{dXx)lJi{dX2) ■ ■ ■ l^(dxn). (5-25) 

Thus the amplitude of the ground state can be written as 

qo{t) = (0o|e-'*^ \4>^)^ j ... j e-^*-ie-^*^^ • • • e-^*^>(dxi)//(dx2) • ■ ■ yi{dxn) 

= (tf ^le-^*-^^ I rf))(#|e-^*^^ I • • . (^J-^le-^*-^" | ^S")). (5-26) 

Therefore the amplitude of the walker at ground state qait) {qk{t) the amplitude of the walker 
at stratum k and time t ) of product graph can be obtained simply by multiplying the cor- 
responding amplitudes of sub-graphs. Also comparing equations (4-17) and (5-26) one can 
determine the spectral distribution ^{dx) and Stieltjes transform G^(x) (consequently all am- 
plitudes of product graph), provided that the product graph possesses quantum distribution, 
too. In the following, we will obtain the required conditions for obtaining QD graphs from the 
product of two given QD graphs. As we will see it is quite possible to obtain QD graphs with 
the same but non-isomorphic spectral distribution (having different adjacency matrices, such 
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that the adjacency matrix of product one may not be tensor product of adjacency matrix of 
multiphcand graphs) from the products of QD graphs. 

Considering two QD graphs, Gi with ujI,uI, ...,0;^; al, al, and G2 with 
ujI,u!2, •■•'^mj Q;2, Q;|, ce^+i, we obtain the following QD graph from their product 

I0o) = 10^)10^), 



102) = ^==_=_(V^i0^)i0^) + 2y;:;Mi0})i0?) + ^Mm), 



|0.nn-l) = 10^)10^), (5-27) 

with 



<^3 = 



CU4 — 



ujIuj2 + 'ioj\uJi + ojIojI 

lOiUjIujI + '^{LOiLoliOi + a;iCJiC<j2) + ^i^l^l 

12 1212 
q;2 = q;2 = q;2, 0:3 = q;2 + 012 = 0:3 = a^, 

= al + al = al + al = al = al, ctfe+i = ka^ (5-28) 

Now, by acting operators A'^,A~ and on the quantum state of the product graph (5-23), 
we get the following required conditions between the QD parameters of multiplicant graphs in 
order to obtain a QD graph from their Cartesian product 

ul + 3ujI = 3ujI + cul, uj\ -\- AuJl = Aijj\ -\- cul, 
3a;] + 12a;? = 4a;3^ + 6a;^, 6a;2^ + 4a;^ = 12a;J + 3a;^, (5-29) 
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Obviously it is difficult to solve the above equations in general. But one can show that these 
equations hold for graphs which are themselves products of some elementary QD graphs with 
only one, two or at most three nonzero QD parameters. 

Thus, in the following, we will consider two important classes of direct product, namely 
class A and class B with some relevant examples for each class. 

Class A. Class A graphs consist of product of an arbitrary numbers of QD graphs with 
only two non- vanishing QD parameters 

cui — a, UJ2 — — ■ • • — 0', ai — 0, a2 — b, — — • ■ ■ — 0, (5-30) 

where for their n-fold direct product we have 

|Vo) = |0o)|0o)---|0o) = |</'o)^" 

= ^ E |0o)|0o)---|0o)|0i), 

V 1 permutation 

IV'2) = ^ E I0o)|0o)---|0o)|0l)|0l), 

V 2 permutation 

|V'n) = |0l)|0l)---|0l) = l0ir" (5-31) 

Application of operators ^4+, A~ and on quantum states of the product graph given above 
yields for its QD parameters: uj'/, = k{n — {k — l))a, for k = 1,2, 3, n, and a'f. = {k — 1)6, 
where the QD property of the product graph imposes the following conditions: | 14 | — 1 > ct^^i 
and = d, where d is an integer. 

Here we give two non-isomorphic graphs out of many prototypes of class A: 
graphl. \ Vk \— uj'iUj'2...,uj'k, {tree) 

graph2. \ V\ |= na, ; | V2 |= 2r? , ; | V3 |= 3!an(n — 2), .... 
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n-cube as an example of class A 

As an example, we obtain n-cube as the direct product of K2 graphs by choosing a — 1, and 
6 = 0. Hence as we will see in the following that its amplitude \ipt) can be written as the 
product of the K2 ones. Since its adjacency matrix can be written as the Cartesian product 
of the K2 graphs therefore 

1 " 

A = - V /2 <8) • • • ® 0-^ (8) • • • ® ^2 (5-32) 

where the ith term in the sum has ax (Pauli matrix) appearing in the ith place in the tensor 
product. Therefore, we have 



exp{-iAt) = /2 ® • • • ® e"'*'"^/" ^ ■ ■ ■ ® h 
1=1 

where S®" is the tensor product of n copies of B. On the other hand for K2 we have 

cui — 1, UJ2 — — • ■ ■ — 0, ai — a2 — ■ ■ • — (5-34) 

and 

= l^^l{S{x-l) + 6{x + l)). (5-35) 
Hence, for the amplitude of observing the walker at 0-th stratum of n-cubc at time t, we have 

qo{t) = iMe-''^ I 0o) = (<)|e-^*-/" | 0«) ■ ■ • (^^-^le^*-/" | cj>^^) 

— cos (t/n) ■ ■ - cos (t/n) = cos" (t/n) (5-36) 

and the other amplitudes can be calculated by the prescription explained above. 

Class B. Class B graphs consist of product of an arbitrary numbers of QD graphs with 
only four non-vanishing QD parameters 

cui — CU2 — a, UJ3 — ■ ■ ■ — 0; ai — 0, a2 — b, — 2b, a4 — — ■ ■ ■ — 0, (5-37) 
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where for their n-fold direct product we have 



l^o) = \(po)\(po)---\<Po) = \<PoY 

1 

permutation 



lV'i) = 7^ E |0o)|0o)---|0o)|0o)|0i 

V 1 permutation 

= , \ (2! E l</'o)l0o)---|0i)|0i))+ E I0o)|0o)---|0o)|02)), 

permutation permutation 



|V'2n) = 102)102) ••• 102) = 102^". (5-38) 

Similarly, application of the operators and on the quantum states of the product 

graph given above yields for its QD parameters: cui^j^^i — {2k-\-l){n — k)a, for k — 0,1, 2, 3, n, 
^2k ~ k{2n — {2k — l))a,for A; = 1, 2, 3, n and ~ nb together with the constraint similar 
to the class A ones plus an extra constraint cui = LO2 coming from its QD property. In the 
following, out of many prototypes of class B, we give the angular momentum graphs as an 
example. 

Example. Angular momentum graphs 

Angular momentum graphs can be obtained from symmetric tensor product of QD graphs of 
class A with a = 2 and 6=1, called vector graphs ( as the angular momentum 1 can be 
obtained from symmetric tensor product of vectors). Again the amplitude \'4)ti can be written 
as the product of vector graphs, since its adjacency matrix can be written as Cartesian product 
of vector graphs 

1 " 

A = - y /7 • • • (8) A' (g) • • • (8) /r (5-39) 

where the ith term in the sum has A (adjacency matrix of vector graphs ) appearing in the 
i-th place in the tensor product. On the other hand, for vector graph we have 

LO\—LOi — 2, LOz — ijJ^ — --- — 0; CKi = 0, a2 — \, 0:3 = 2, q;4 = q;4 = • • • = 

(5-40) 
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and 

(5-41) 

Hence, for I — n angular momentum graph with the following QD parameters 

uji^2n, a;i = 2(2n- 1), • • -,0;^ = A;(2n - A; + 1), • • •; ctfe = A; - 1, (5-42) 
the amplitude of observing the walker at 0-th stratum of at time t becomes 

= le-^*/"(2 + 3cos(V5i/n) + iV^sin(V5Vn)) • • • \e-'*'''{2 + 3cos(V5i/n) + iVh sm.{Vbt / n)) 
5 5 

= (ie-**/'^(2 + 3 cos(\/5t/n) + i\/5 sin(\/5t/n)))", (5-43) 
5 

and the other amplitudes can be calculated by the prescription explained above. 

6 Some examples of QD graphs 

In this section we provide some examples of finite and infinite QD graphs and using the spectral 
distribution we calculate the relevant amplitudes of continuous time quantum walks on theses 
graphs. As far as the authors know, theses examples exhaust all known graphs which have 
been studied by researcher in quantum walk until now. First we begin with finite QD graphs: 

6.1 Complete graph Kn 

A complete graph with n vertices (denoted Kn) is a graph with n vertices in which each 
vertex is connected to the others (with one edge between each pair of vertices) . Therefore, it 
is a two-stratum QD graph with one non-vanishing QD parameter uji = n — 1 and spectral 
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distribution 

T — n — I— 2 n — 1 1 

G,{x) = ^- -f — , 6{x+l) + -6{x-{n-l)), (6-44) 

— [n — 2)x — n + 1 n n 

that yields the following amplitudes at time t of the quantum walker 

JR 

^ !i^e^ + -e-'' = -(e-^* + (n - l)e^) (6-45) 

qJt) = (0i|e-^^*/("-^) I 0o) = / ^—e-'^''^''-^^nidx) 

jRn — 1 

^/^^(e-*-e^*/("-i)), (6-46) 



where the results thus obtained are in agreement with those of Ref.[15]. 
6.2 Cycle graph Cn 

A cycle graph or cycle is a graph that consists of some number of vertices connected in a closed 
chain. The cycle graph with n vertices is denoted by C„, and quantum walk on them turns 
out to be different for odd and even n, hence we treat them separately. 

Odd n. For odd n = 2m + 1, C2m+i graph consists of m -|- 1 starata, and its QD parameters 
are 

a;i = 2, tc;2 = = • • • <^m = 1; ai = q;2 = • • • = 0, a^+i = 1. (6-47) 
It is straightforward to show that it has the following spectral distribution 

-1 r) 2m 9/-7r 

= —^S{x - 2) + -^y^5{x - 2cos(-^^)). (6-48) 
Therefore, the amplitude probability for walker at time t and 0th stratum is 

JR 



1 2 



m 



-it _L ^ ^ g-it cos 2i7r/(2m+l)_ (6-49) 



2m -M 2m -Fl 
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Using Eq.(4-17), one can calculate the other amplitudes, where the results obtained are in 
agreement with those of Ref.[15]. 

Even n. For even n — 2m, graph consists of m + 1 starata, and its QD parameters are 
a;i = 2, a;2 = a;3 = • • -(^m-i = 1; '^m = 2; ai = 02 = • • • = 0. (6-50) 
One can straightforwardly show that it has the following spectral distribution 

-1 o 2m— 1 Ol-TT 

W.-2, + .(. + 2))+ E .(.-2cos(-)). (6-51) 

t=l,l^m 

Similarly, the amplitude probability for walker at time t and — th stratum is 



R 

1 1 ™ 

^ _Lcost + — y e-**^°^2'"/(2'"). (6-52) 

Again one can calculate the other amplitudes by using Eq.(4-17), where the results are in 
agreement with those of Ref . [15] . 



In the limit of large n, the roots xi = 2cos(|^) reduce to xi = 2cos(7ra;) with x — 
lim^^;_^oo ^ and both spectral distributions given in (6-48) and (6-51), irrespective of whether 
n is even or odd, reduce to continuous spectral distribution ^{x)dx with ^{x) 

1 r 



H{x) = - / dyS{x - 2 cos(y)) 

TT Jo 

n Jo 2 sin(y) 

1 ^ S{y — aYccos{x/2)) 



TT Jo 2 sin(y) 

1 1 



-2 < a; < 2, (6-53) 



TT V4 - x2 

which is the same as the continuous spectral distribution of infinite line graph Z. 
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6.3 Tchebichef graphs 

By choosing Tchebichef polynomials of first kind (second kind )with scaling factor ^ as or- 
thogonal polynomials appearing in recurrence relation (4-11), i.e., Qn{x) — 2('"~^^"+^T„(a;/2"^) 
(2('"~^)"C/„(x/2"*)), one can obtain a class of finite and infinite QD graphs of Tchebichef type, 
with QD parameters oui = 22(™-i)+i, = 22('"-i), A; = 2, 3, and ak = 0, k = 1, 2, 3, ... 
{ujk — 2'^^"^~^\ k — 1,2,...; ak — 0, A; = 1, 2, ...), such that its Stieltjes transform of spectral 
distribution becomes GJx) = ,f"^rA ). 

Therefore its spectral distribution can be written as = ^Z)/=o^(^ ~ 2"^ cos ^^^2n^^ ) 
E;=o sin2( J^)5(x-2- cos ^)), where 2- cos Z = 1, 2, ... (2- cos(;ff^), Z = 1, 2, • • •) 

are roots of the Tchebishef polynomial, T„(x/2™) (f/„+i(a;/2'")). Now, using Eq.(4-17) the am- 
plitude of observing walker at different stratum at time i can be calculated as 

1 -r^ _j2'"tcos(i2i+ilZL) 



for the first kind and 



q,{t) = ^ Vsin2(-^)e-^2'"*^°«(^), 



q^^t) = ^ Vsin2(-^) (A: + l)/7r _,,^,,,,(_^) ^ ^ _^ 
for the second kind. 

Finally, in the limit of large n, the absolutely continuous part of spectral distribution ij,{x) 
reads as Ijl{x) — (^^2^"* — x^) such that, —2"* < x < 2'"(its derivation is similar to 

that of cycle graph C„) and the amplitude Pk{t) given in (6-54) and (6-55) reduce to 

qo{t) = Jo(2— 4), 

qk{t) = , / 2*^(— i)+irfc(a;/2™)e-^^*/V((^a;) 
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= \/2i'=Jfc(2™-^t), k>l, (6-56) 

for the first kind and 

qk{t) = ^ [ 2^'^-^^^Uk{x/r')e-'^'/^n{dx) 

27r ^1 - {x/2^'^y 

TT Jo 

= i\Jk{2"'t) + Jk+2{2rt)), A; = 0, 1, 2, (6-57) 

for the second kind, where Jkif) are the Bessel polynomials. In deriving of (6-56) and (6-57), 
we have used the integral representation of Bessel function after making the change of variable 
X — 2"*cos'0- Also, one can calculate the amplitudes of infinite Tchebishef graphs by using 
generating function of Tchebishef and formula Eq.(4-r7), where the obtained results are in 
agreement with those given in (6-56) and (6-57). As we will see in the following, some of the 
known finite and infinite graphs can be obtained from Tchebishef graphs by appropriate choice 
of m and the polynomials as: 

A. For m — 1 and Tchebishef polynomial of the second kind we obtain finite path graph 
Pn — {0, 1, 2, ...}, where it is a n- vertex graph with n — 1 edges all on a single open path. Its 
QD parameters and the amplitude of observing walker at different stratum at time t are 

= a;2 = a;3 = • • • = 1; = 

q^it) = ^ Vsin2( J!L-)e-*-^(^), 

qk{t) = ^ sin(i7r/(n + 2)) sin(^^Z7r)e-^*^°^(''^/("+2)), k>l, (6-58) 
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where the results thus obtained are in agreement with those of Ref.[26]. Also, in the limit of 
large n, qk{t) reduces to 

Qkit) = i'iJkit) + Jk+2(t)), k^O,l, .... (6-59) 

B. For m = I and Tchebishef polynomial of the second kind we obtain a sequence of 
graphs Gn (for more details see Ref. [7]). The number of vertices in Gn is 2""*"^ + 2" — 2. In 
Fig.l we show G4. In general, G„ consists of two balanced binary trees of depth n with the 
2", nth-level vertices of the two trees pairwise identified. For the quantum walk on G„, one 
starts the waker at the root of a tree and wants to calculate the probability of the presence of 
the walker at the other vertices as function of time. One can show that G„ is QD graph with 
(2n -I- 1) strata, where stratum j consists of 2^ vertices for j — 1, 2, n -|- 1 and 2^^"+^"^^ for 
j = n + 1, ...,2n + 1. Therefore, its QD parameters and the amplitude of observing walker at 
different stratum at time t are 

a;i = a;2 = a;3 = • • • = 2; cci = 0:2 = • • • = 0, 

2 ^ . 2(_^^^-i2V2tcos{^) 



= :-^Esin'(— -r)e 



9 t -I- 1 r- 

q^{t) = ^— J2 sin(Z7r/(n + 2)) sin(^^/7r)e-''^*'^°^(''^/("+2))^ ^ > 1_ (6.60) 

In the limit of large n, Pfe(i) reduce to 

q^{t) = ^^=(7^(2^2^) + Jfe+2(2V2i)), A; = 0, 1, .... (6-61) 

. C. For m — \ and Tchebishef polynomial of the first kind we obtain infinite line graphs Z 
where its QD parameters and the amplitude of observing walker at different stratum at time 
t are 

a;i = 2, a;2 = a;3 = • • • = 1; = 0, 

?o(i) = Jo(i), 
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quit) = yJ^t'^Jkit), k>l. 
The results thus obtained are in agreement with those of Ref . [27] . 
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(6-62) 



6.4 The finite Hermit graph 

Let G consist of n vertices such that its QD parameters, cardinahty of stratums and Stieltjes 
transform are 

cui — n, UJ2 — n — 1, CU3 — n — 2, ■ ■ ■ ,u!n — i, ai — a2 — ■ ■ ■ — 0, 



Vi\^n, \V2\^n(n-l),---,\Vn\^n\, 
where, Hn{x) is Hermite polynomials. Then it will have the following spectral distribution 



GM = (6-63) 



pi^^JZ^i^-^i)^ (6-64) 

n-\- I ^ 

where xi, I — 1,2, ... are the roots of Hermite polynomial Hn{x). Again, using Eq.(4-17) or 
the recurrence relation (4-18) one can calculate all amplitudes, where in the following we just 
quote the results for n = 3. Its QD parameters and spectral distribution are 



a;i = 3, a;2 = 2, cua = 1; cti = q;2 = ■ ■ ■ = 0, 



and the amplitude as 



qo{t) = ^(cos(\/3 + V6)t + cos(V3- \/6)t). 



qiit) = -^((V3 + V6)sin(V3 + V6)t + (V3- V6) sin(V3 - V6)t), 

q2{t) = -^{cos{\/ 3 + V6)t - cos{'\/3-V6)t), 
2v6 



q^{t) = -^(^3 + V6(V3-V2)sin(\/3 + V6)i+\/3- V6(V3 + V2)sin(\/3- V6)i). (6-66) 
6v2 
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6.5 The infinite Hermit e graph 

By choosing Hermite polynomial Hn{x), as orthogonal polynomial appearing in recurrence 
relation (4-11), one can obtain a class of infinite QD graphs of Hermite type with QD param- 
eters cuk — k, and ak — 0. Here we give two non- isomorphic graphs out of many prototypes of 
infinite Hermite graph: 

Graphl. \ Vk \— <^i<^2 ■ ■ - ^^k — k\{tree)i.e., K^{y) — 1 

Graph2. \ h 1, 1 h 2(«:_(?/) = 1), | h 3!(«_(y) = 1), | h 3!(«-(y) = 2), | h 
30(«:_(y) = 1), I 1^6 1= 6 X 30(«:_(y) = 1), .... 

Here the amplitude qk{t) can be calculated rather trivially by using generating function of 
Hermite polynomials as follows: 

^ V 2/c!7r Jr ^ k\ J-oo 

= e-y'/^ r e-^/2(^'+2"M)dx = e-^*(*+2^^). (6-67) 

J —oo 

Now, by comparing the terms of the expansion appearing on the left hand side of (6-67) with 
those of expansion on the right hand side we obtain 

quit) = ^^e"*'/'' (6-68) 

where the result thus obtained is exactly the same as that of WKB approximation, i.e, the 
method of stationary phase gives the exact result in infinite Hermite case. Then, the probability 
distribution | qk{t) '^^^ ^^ commonly called the Poisson distribution. 

6.6 The infinite Laguerre graph 

We take the Q„(x) as a Laguerre polynomial, such that they satisfy the orthogonal polynomials 
appearing in the recurrence relation (4-11), i.e, Qn{x) = n\{—lYa^L^\'^). Then, one can 
obtain a class of infinite QD graphs of Laguerre type , with QD parameters ujk — a^k{k + 7), 
and choose b — —0(1 + 7), ai = and a^+i = 2ka ior k > 1. Here we give two non- isomorphic 
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graphs out of many prototypes of infinite Hermite graph: 

Graphl. \ \— UJ1UJ2 ■ ■ ■ u!k{tree) i.e.,K^{y) — 1 

Graph2. \ 14 |= k\a'^{l + 7)(2 + 7) ■ ■ ■ (A; + 7), = a). 

By using Eq.(4-17), we obtain generating functions for amphtude probabihty qk{t) 



k=0 



- - r e- \~'^^A ^re-i^^dx = , f^^,,, . (6-69) 

a Jo (l-y)T+i^ a ' (l + ia(l-y)i)T+i ^ ' 

It is straightforward to show that Pk{t) has the following form 



Formula (6-70) implies that for distant time the amplitude becomes proportional to the inverse 



of time t, i.e, we have qk{t) ~ y Ity!^' |y?t, which is the same as the one obtained by WKB 
approximation. 

7 Asymptotic behavior of continuous-time quantum walk 
in the large time limit 

7.1 Average moments of number of visiting strata 

Studying the large time behavior of moments and variance as function of t for classical random 
walk and quantum walk is important, since common properties for describing the behavior of 
a walker are the exponents u, characterizing the scaling with time t of the mean square end-to- 



end distance (i?^) = y {k"^) — (/c)^ ~ t''(the exponent u determine universality class). The walk 
on infinite line, | qk{t) P has the form of a binomial distribution, with a width which spreads 
like i^/^(i.e, u — 1/2), therefore the variance grows linearly with time. But the variance in the 
quantum walk on infinite line, by contrast, grows quadratically with time, and the distribution 
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I qk{t) p has a complicated, oscillatory form[28]. Now, by using the probability amplitudes of 
Eq.(6-68) and Eq.(6-70) we evaluate the average of different moments of stratum number, i.e, 
{k'^) (g-th moment, q = 1,2,3, ... ), for quantum walk on infinite Hermite, Laguerre and line 
graphs as follows: 

A. Using the generating function Hermite polynomials, the average g-moment of infinite Her- 
mite graphs can be written as 

Above formula implies that (k) = and (fc^) = + t'^, therefore for it standard deviation 
a{t) we get a{t) = ,J{k^) - {k)^ = t. Also, it is straightforward to see that for large time, the 
g-moment has the following power dependence upon time 

(A;«) {ty. (7-72) 

B. Similarly for infinite Laguerre graphs we have 

/;,.^ f ;,. (7 + fc)! jatr _ 1 .d 1 

^ '~h '^'t! (l + (at)2)fc+7+i - (l + a2t2)7+i^^rf^^ ^(1-^)7+1^1.=^' 

where it implies that 

(A;) = (7 + l)aV, (A;^) = (7 + l)aV + (7 + 1)(7 + 2){aH''f 

lim {a{t)) = ^^ lim (A;*) (a^t^)^. (7-74) 

t— >oo t— >oo 

C. In infinite line graph, where the walker starts at origin, we calculate the average moments 
of visiting number of sites (therefore k denotes the sites of graph where k — Q correspond to 
its origin), again using the well known the generating Bessel function we can show that 

00 

{k) = E I Jk{t) |'= 0, 

A;=— 00 

00 

k=—oo 
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for q =even 



otherwise. 



a{t) t. 



(7-75) 



In the above examples the the average number of visiting strata (sites for infinite hne), the 
quantity, ait) is the same as end-to-end distance. The end-to-end distance for quantum walker 
on infinite Laguerre graphs varies quadratically with time, i.e, (i?^) ~ t'^, while the end-to- 
end distance for walker on infinite Hermite graphs is (i?^) ~ t. Thus, the quantum walk 
on infinite line and Hermite graphs have the same critical exponent, which is different from 
infinite Laguerre one. Therefore, quantum walk on infinite line and Hermite graphs belong to 
the same universality class, which is different from the universality class of quantum walk on 
infinite Laguerre graph. 

7.2 Asymptotic behavior of the probabihty distribution in the large 
time limit 

In general we can not have an analytic expression for the amplitudes of continuous time quan- 
tum walk on most of graphs , i.e, the integral appearing in the Eq.(4-17) is difficult to evaluate, 
but one can approximate it for large time by using the method of stationary phase explained 
above in section 4. In order to see the importance of spectral distribution method, we obtain 
the asymptotic behavior of some finite and infinite graphs in remaining part of this section. 

7.2.1 The infinite graphs: 
7.2.1.1 Star lattice 

Star lattice is an A^-fold star power G*^ , where G = (V, E) is the half line of integer, i.e., V — 
0, 1, 2, ... and i ~ j' if and only ii \ i — j \— 1, so its QD parameters are cui — N,u!k — i,Oik — 
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for A; = 2, 3, and as it is shown in Ref.[19] it as the following spectral distribution 

Therefore, its ground stratum amplitude qo{t) can be written as 



Now, its asymptotic form can be obtained by the method of stationary phase after making the 
change of variable x — 2cos(^), where we get 

, , 4iVr(3/2) 1 371, 

lox N ^ 2, while, for A?" = 2 the star lattices reduces to infinite hne graphs with spectral 
distribution = :^;yj=f, — 2 < a; < 2 and asymptotic ground amplitude qoit) = Jo{2t) pa 

7.2.1.2 Two-dimensional comb lattice 

As another example we consider two-dimensional comb lattice of Ref. [29], with the fol- 
lowing spectral distribution 

fjL{x)dx ^ 1 ^ -2v^<x<2v^ (7-79) 

TT V8 — 



and ground stratum amplitude 

q,(t) = / e-'''/'p{x)dx = / L^^J^ 
J-2V2 J -2V2 n 8 — x'^ 

1 r 



TT Jo 



Again its asymptotic form can be obtained by the method of stationary phase after making 
the change of variable x — 2-\/2cos(^), where we get 



qo{t) - y^cos(t/y2-7r/4). (7-81) 
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7.2.2 The finite graphs: In order to obtain the asymptotic form of finite graphs at large 
times , we need to use Gauss quandrature formula(4-21) to approximate finite graphs with the 
infinite ones and then using the method of stationary phase we can obtain their asymptotic 
form. In the example under discussion asymptotic behavior of quantum walk on line and 
generalized to asymptotic behavior quantum walk on finite graphs. For line qo{t) is Jo{t), then 
for large time t we have qo(t) ^ ^J^cos{t — f )• If we use (4-21) and calculate go(^) on finite 
graph line we have qo{t) = i ELo e-^*^°^(^'^). Now, we calculate numerically the difference 
of amplitudes of infinite line with finite graphs 7r{n, t), for large time t 



The results are depicted in Fig. 2 and Fig.3, where 7r(n, t) is in the limit n from 500 to 600 and 
t in the limited 1000, is limited to zero. Then, to study the behavior of asymptotic quantum 
walk on finite graph line, we study the behavior of asymptotic quantum walk on infinite graph 
line. Therefore to study the behavior of the asymptotic quantum walk on finite graphs, we 
can use arithmetic, approximate it with infinite graph, and by using the method of stationary 
phase, study the behavior of asymptotic quantum walk. 

8 Conclusion 

Using the spectral distribution associated with the adjacency matrix, a new formalism for 
investigation of continuous-time quantum walk on some graphs is developed, where the or- 
thogonal polynomial together with their recursion relations play an important role. As as by 
product, by using the Gaus quadrature formula one can approximate infinite graph with finite 
one and vice versa, which leads to the derivation of large time asymptotic form of continuous- 
time quantum walk amplitudes of both finite and infinite graph, simply by using the method 



Tr{n,t) =1 Jo{t) ^ e 



-j^ n—l 



itcos(2|±l7r) 



k=0 




(7-82) 
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of stationary phase. Even though the powerful method of spectral distribution seems to work 
for some restricted kinds of graphs introduced in Ref. [19], but it is possible to generalize it to 
work for most of graphs, particularly for the graphs with strata dependent {y) and regular 
distant ones, which is under investigation. 



In this appendix we prove the following lemma in connection with the equality of continuous- 
time quantum walk amplitudes on the vertices belonging to the same stratum. 

Lemma 1. Let qik{t) denote the amplitude of observing the continuous-time quantum 
walker at vertex i e Vfc at time t. Then for a class of QD graphs, the amplitude qik is the same 
for all vertices of stratification k, for all t. 
Proof. 

Let as take the Fourier transform of unit vectors |A;, i) of the stratum Vk'. 



where u!\Vk\ = s'^''' • Now, acting the adjacency matrix A on it according to the formula (3-4), 
we obtain 
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1 



\Vk\-l 




(A-i) 



A\(/)k,o) = V^fc+i|0fc+i,o) + oik+i\(t>k,o) + V^l^fc-i,o)- 



(A-ii) 



Therefore, {(f)m,i\A 



0fe,o) — 0, for I 0, hence we have 




1 



Qk{t) 



Pokjt) +Plk{t) +P2k{t) H ^P(|Vfc|)fc - 1)(^) 

\yk\ 



for,l = 0, 



(A-iii) 
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or 



qok{t) 
qik{t) 

(l2k(t) 
\ (l{\Vk\-l)k{t) j 



^ qk{t) ^ 








V 







(A-iv) 



where | | x | | the discrete Fourier transformation matrix(DFT) F is defined as 



1 


1 


1 


... 1 


1 


cu 




... J^"l-1 


1 






^2(|y„|-i) 



^ 1 a;2(|y„|-i) 



Inverting F in Eq.(A-iv) we obtain gofc(^) = qik{t) = • • • = q{\Vk\-i)k{t) = 



(A-v) 



Appendix B 

list of orthogonal polynomials connected with some finite and infinite graphs 

In this appendix we give a list of orthogonal polynomials (in addition to those mentioned 
above in section 6), where after appropriate one-dimensional, affine transformation(shift and 
recalling), they can be reduced to monic orthogonal polynomials with recursion relation(4-ll), 
such that the corresponding parameters a;„ and can be connected with some finite or infinite 
graphs of QD type. 

l.Chcirlier polynomials 
Qr,{x) = a"C('^)(^), fi{x) = j:^^,sZ^6{x-j), interval - (0,+oo), b = 
—ad, ujn = ch^nd, ot^+i = na. 
Graphl.j Vn \— ujxloi ■ • •a;„(tree), 
Graph2.| Vn |= n\o?d^ and = a. 
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2. Meixner polynomials of the second kind 

Qn{x) = a"M„(2^), i^{x) = n{x;5,r)) = {r{r)/2))-^ \ r{^) p e-^*^'^"'^'^), interval = 
(— oo,+cxo), b — —ari6, Un = a^n{n + r] + 1){6^ + 1), an+i = ^.anS. 
Graphl. | Ki |= UJ1UJ2 • • •(^n(tree), 

Graph2. | K |= n!a^(2 + ri){3 + rj) ■ ■ ■ (n + 1 + rj) and K-(y) = a. 

3. Orthogonal polynomials related to Jacobi Elliptic functions 
3.1 

Qn(x) — a'^An(^): n{x) — 7r2(2n)!(2^+i)!fc2n '^t-qi'^^+'^) interval — 

(0,+oo), b = -a{l + k'^), ujn = '^a^n^{^n^ -l)k'^, On+i = 4an(n + 1)(1 + A;^). 
Graphl. | |= UJ1UJ2 • ■ •a;„(tree), 

Graph2. | K |= A''{n\fa'^{l x 3 x 15 x ... x (471^ - 1))A;2" and = a. 

3.2 

g„(x) = a"S„(^), = Ejtg^ 2,r^(2n+i)!f2l+2)!fc^" . ■'^(^^-j)^ mtcrmZ = (1, +oo) , 6 = 

-4a(l + A;2), a;„ = aMn(n + l)(2n + 1)^/^^ q;„+i = 4an(n + 2)(1 + A;^). 
Graphl. | Ki |= OJ1OJ2 • ■ -(^nltree), 

Graph2. | K h 4"n!(n + l)\a^{2,'^ x 5^ x ... x (2n + l)^)^;^" and = a. 

3.3 

= Cn{x), fi{x) = E/=^oo fcj^(n)!2fc2[n/2j i+q(tj+i) Kx-j), interval = (-00, +00), b = 
-4a(l + A;2), a;2n+i = (2n + l)^ u;2n = (2n)2A;2, = 0. 

Graphl. | |= '^i'^2 • • -o^nltree). 
3.4 

Q„(a;) = Dn{x), ^{x) = Ej=-^ K{n)Pk^i(^+i)/^] i+g2, '^(a^-i), interval = (-00, +00), u;2n+i = 
(2n+l)2/c2, cu2n = (2n)2, a„ = 0. 
Graphl. | 14 |= UJ1CU2 ■ • •a;n(tree). 

Such that, j = 0, ±1, ±2, ■ ■ ■ , K = /o°°(l - ^ sm\(j)))-^/^d(j) 

and is a certain constant related to k which appears in the theory of theta function. 
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3.5 Carlitz polynomials 
3.5.1 

Q^(x) = a"F„(2^), b = -a, = Ao?n^{2n - 1)^^, a^+i = 4an(n(l + P) + 1). 
Graphl. | Ki |= UJ1UJ2 • • •(^n(tree), 

Graph2. | K h 4"a2n!2(l x 3^ x 5^ x ... x (2n - 1)2)A;2" and = a. 

3.5.2 

g„(a;) = a'^G'„(3^), b^-ak'^, ujn^^a?n'^{2n-lfk'^, q;„+i = 4an(n(l + A;^) + A;^). 
Graphl. I \^ |= UJ1UJ2 ■ ■ •Ci;„(tree), 

Graph2.| I4 |= 4Vn!2(l x 3^ x 5^ x ... x (2n - l)^)^^" and H-{y) = a. 
3.5.3 

Qn(a;) = a"G*(^), 6=-a(4 + A;2), = 4aV(2n+ 1)2A;2, q;„+i = 4an(n+ 1)(1 + A;^). 
Graphl. | |= 1^11^2 " " •<^n(tree), 

Graph2. | K h ^""^^{1 x 3^ x 5^ x ... x (2n - 1)2)A;2" and = a. 

The explicit weight function for the subsection (3.5) can be obtained routinely from the weight 

function of the subsections (3.3) and (3.4). 
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Figure Captions 
Figure- 1: Figure G4. 

Figure-2: Shows 7r(n, t) for the infinite line graphs with the n-vertices finite graph (obtain 
from the infinite on via Gauss quadrature rule) as a function n att — 1000, where the difference 
is almost negligible for n > 518. 

Figure-3: FigurcS show 7i{n,t) for the infinite line graph with the n-vertices finite graph 
as a function of time and number of vertices n, where the difference is almost disappear for 
t > 1000 and n > 518. 



